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if (env=getenv("_OGL_GAMMA_FILTER")){ ■'■ 
/* Gamma filtered*/ 1 0P° 

float gamma; f 

float sum; 

int factor; 

int width; 

int size; 

gamma=1.0; 
gamma=atof(env); 

width=8; 

if (env=getenv("_OGLGAMMA_TABLEWIDTH")) 

width=atoi(env); 
factor=(int)pow(2.0, (double) (8.0-width)); 

sum=256.0; 

if (env=getenv("_OGI GAMMA_SUM")) 

sum=atof(env); 

size=256; 

if (env=getenv("_OGl_GAMMA_TABLESIZE"))| 
size=atoi(env); 
switch (size) { 
case 16: 

for (i=0; i<256; i+=16) j 

a=sum ■* pow((double)(i/256.0), (double)(1.0/gamma)); 
b=sum * pow((double)((i+1)/256.0), (double)(1. 0/gamma)); 
c=sum * pow((double)((i+2)/256.0), (double)(1. 0/gamma)); 
d=sum * pow((double)((i+3)/256.0), (double)(1. 0/gamma)); 
e=sum * pow((double)((i+4)/256.0), (double)(1. 0/gamma)); 
f=sum * pow((double)((i+5)/256.0), (double)(1 .0/gamma)); 
g=sum * pow((double)((i+6)/256.0), (double)(1.0/gamma)); 
h=sum * pow((double)((i+7)/256.0), (double)(l. 0/gamma)); 
i=sum * pow((double)((i+8)/256.0), (double)(l. 0/gamma)); 
j=sum * pow((double)((i+9)/256.0), (double)(l .0/gamma)); 
k=sum * pow((double)((i+10)/256.0), (double)(1. 0/gamma)); 
l=sum * pow((double)((i+11)/256.0), (double)(1 .0/gamma)); 
m=sum * pow((double)((i+12)/256.0), (double)(1. 0/gamma)) 
n=sum * pow((double)((i+13)/256.0), (double)(1. 0/gamma)); 
o=sum * pow((double)((i+14)/256.0), (double)(1. 0/gamma)); 
p=sum * pow((double)((i+15)/256.0), (double)(1. 0/gamma)); 
AARIterTable[i]=AARIterTable[i+1]= 
MRIterTable[i+2]=AARIterTable[i+3]= 
AARIterTable[i+4]=AARIterTable[i+5]= 
AARIterTable[i+6]=AARIterTable[i+7]= 
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AAniterTable[i+8]=AAniterTable[i+9]= 6 / 7 ; -V 

MRIterTable[i+10]=AARIterTable[i+1l]= 

MRIterTable[i+12]=AARIterTable[i+13]= 

MRIterTable[i+14]=AARIterTable[i+15]= FIG 1 OB % " ' 

(int)((((a+b+c+d+e+f+ 
(int)((((a+b+c+d+e+f+ 
g+h+ii+j+k+m+ 

n+o+p)/(16.0*factor))+0.5)*factor); 

! 

break; 
case 32; 
for (i=0; i<256; i+=8) { 
a=sum * pow((double)(i/256.0), (double)(1 .0/gamma)); 
b=sum * pow((double)((i+1)/256.0), (double)(1 .0/gamma)); 
c=sum * pow((double)((i+2)/256.0), (double)(1.0/gamma)); 
d=sum * pow((double)((i+3)/256.0), (double)(1.0/gamma)); 
e=sum * pow((double)((i+4)/256.0), (double)(1.0/gamma)); 
f=sum * pow((double)((i+5)/256.0), (double)(1.0/gamma)); 
g=sum * pow((double)((i+6)/256.0), (double)(1.0/gamma)); 
h=sum * pow((double)((i+7)/256.0), (double)(l.O/gamma)); 
AARIterTable[i]=AAniterTable[i+l]=AAniterTable[i+2]= 
AAniterTable[i+3]=AAniterTable[i+4]=AAFilterTable[i+5]= 

MRIterTable[i+6]=MRIterTable[i+7]=(int)((((a+b+c+d+e+f+g+h)/(8.0*factor))+0.5)*factor); 

! 

break; 
case 64; 
for (i=0; i<256; i+=4) { 
a=sum * pow((double)(i/256.0), (double)(l.O/gamma)); 
b=sum * pow((double)((i+1)/256.0), (double)(1.0/gamma)); 
c=sum * pow((double)((i+2)/256.0), (double)(1.0/gamma)); 
d=sum * pow((double)((i+3)/256.0), (double)(1 .0/gamma)); 
AARIterTable[i]=AAniterTable[i+ 1 ]= 
AARIterTable[i+2]=AARIterTable[i+3]= 
(int)((((a+b+c+d)/4.0*factor))+0.5*factor); 

I 

break; 
case 128; 
for (i=0; i<256; i+=2) { 
a=sum * pow((double)(i/256.0), (doubie)(1 .0/gamma)); 
b=sum * pow((double)((i+1)/256.0), (double)(1 .0/gamma)); 
AARIterTable[i]=AARIterTable[i+l]= 
(int)((((a+b)/2.0*factor))+0.5*factor); 

i 

break; 
case 256; 

for (i=0; i<256; i++) { 

AARIterTable[i]= 

(int)((((sum * pow((double)(i/256.0), (double)(1.0/gamma)))/factor)+0.5)*factor); 

I 

break; 

\ 
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Assumptions: Floating point coverages are defined in the normalized 0.0 
to 1.0 range in which 0.0 corresponds to no coverage and 1.0 corresponds 
to full coverage. Fixed point coverages are defined in the range 0 to size 

float * GenFloatingPtGammaTable(int size, 

float gamma) 

i 

int i; 

float * table; 

if (table=malloc(sizeof(float)*size)) { 
for (i=0; i<size; { 
table[i]=(float)pow((double)i/(size-1 ), (double)(1 .0/gamma)); 

I 

i 

return (table); 

i 

int * GenFixedPtGammaTable(int size, 

float gamma) 

i 

int i; 

int *table; 
float val; 

if (table=malloc(sizeof(int)*size)) j 
for(i=0; ksize; i++) { 
val=(float)pow((double)i/(size-1), (double)(1.0/gamma)); 
table[i]=(int)((size-1)*val+0.5); 

! 

! 

return (table); 

i 



